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Abstract 

We study in this paper the possible occurrence of orbital magnetim for two- 
dimensional electrons confined by a harmonic potential in various regimes of temper- 
ature and magnetic field. Standard coherent state families are used for calculating 
symbols of various involved observables like thermodynamical potential, magnetic 
moment, or spatial distribution of current. Their expressions are given in a closed 
form and the resulting Berezin-Lieb inequalities provide a straightforward way to 
study magnetism in various limit regimes. In particular, we predict a paramagnetic 
behaviour in the thermodynamical limit as well as in the quasiclassical limit under 
a weak field. Eventually, we obtain an exact expression for the magnetic moment 
which yields a full description of the phase diagram of the magnetization. 

PACS: 75.20.-g: diagmagnetism, paramagnetism, superparamagnetism, 75.30.Kz: 
magnetic phase boundary, 71.10.Ca: Electron gas, Fermi gas, 51.60: magnetic phenomena 
in gases 

I Introduction 

In a recent paper, Ishikawa and Fukuyama []]] describe the possible orbital magnetism 
for two-dimensional electrons confined by a harmonic potential in various regimes of tem- 
perature and magnetic field. They afford a quite large complement of information in 
regard to the previous paper p| devoted to the same subject. According to the range 
of values assumed by the relative ratios between the three characteristic energy scales 
present in the model, namely the thermodynamical unit ksT, the magnetic quantum 
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huj c , and the harmonic quantum hujo, they explain the existence of the different mag- 
netic regimes. As a matter of fact, they distinguish between the "Mesoscopic Fluc- 
tuation" regime (ksT ~ h(^Juj 2 + Aujq — u c )/2), the "Landau Diamagnetism" regime 

(k B T ~ h(^u 2 c + Aul + uu c )/2) , and the "de Haas -van Alphen" regime. Their studies rest 
upon the derivation of an approximate closed formula for the thermodynamical potential Q 
and from which they are able to get the magnetic moment. The isotropic two-dimensional 
harmonic potential offers us a complete analytical treatment of the problem whereas the 
orbital magnetism is expected not to depend strongly on the shape of the confinement. 
The crucial point in our derivation is that the Fermi-Dirac function is a fixed point of 
the Fourier transform. Exact series expansions ensue by simple application of the residue 
theorem as shown in the appendix. Our results for the thermodynamical potential, the 
orbital magnetic moment, and the spatial distibution of currents are easily tractable and 
analyzable in comparison with those given in Ref.0. 

One can show that the Hamiltonian under a perpendicular magnetic field can be 
decomposed into the sum of two harmonic oscillators with frequencies uj + and respec- 
tively. This fact makes the present model completely integrable and also makes the use 
of the coherent states (CS)|| possible. These coherent states are defined as tensor prod- 
ucts of two standard one-dimensional coherent states and read as: {| Z4, z g ) — \ Za ) ® \ 
z g ), z d , Zg G C}. This family of states provides us with a resolution of the unity and in 
this way we are allowed to use the so-called symbol calculus a la Berezin-Lieb-Perelomov 
f\ Consequently, the Berezin-Lieb inequalities give us the opportunity to explore the be- 
havior of the magnetization of our model in extreme conditions like the thermodynamical 
limit or the quasi-classical limit. 

In the next section, we shall give a short review about the physical model. The eigen- 
values and eigenvectors of the Hamiltonian will be obtained by algebraic and analytical 
methods respectively. In particular, we shall explain the symmetries of the problem and 
their relevance in various approximations of the model. In Section III, the construction of 
the coherent states is presented and the main CS algebraic and functional properties are 
briefly recalled there. Section IV is devoted to the thermodynamical potential Q. First we 
apply the Berezin-Lieb inequalities and study the magnetic moment in the different limit 
regimes. We then give the exact expression of Q. The core of the paper lies in Section V, 
in which we establish exact expressions for the magnetic moment and the average number 
of electrons and discuss their behavior in different temperature and magnetic field regions. 
We show in Section VI how coherent symbol calculus yields a Fourier integral expression 
for the radial distribution of the current. This integral can be also given as an expansion 
series with the aid of the residue theorem. In the conclusion we shall give some remarks 
and comments on the possible extension of our approach to the other systems of physical 
interest. 



a - ) Sec [Q and Licb in 0. For a recent up-to-dated review on CS, specially on their group representation 
aspects we also refer to Notice that here the operator calculus in Hilbcrt space is replaced by the 
functional calculus involving CS parameters. 
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II Hamiltonian and Symmetries 



The Hamiltonian for two-dimensional spinless electrons confined by an isotropic harmonic 
potential and submitted to a constant perpendicular magnetic field is written as 

w =i (p+ ; A)2+ H R2 - <"> 

where Coulomb interactions are neglected. This model is called the Fock-Darwin Hamil- 
tonian in the litterature 0. The radius of the system, R m , is classically defined as 

l -mulR 2 m = //, (2.2) 

where \x is the chemical potential, and we shall give in Section IV a full quantum statistical 
mechanics interpretation of this relation. We shall work with the symmetric gauge A = 
|H x R. We first solve the eigenvalue and eigenvector problems by using the underlying 
Weyl-Heisenberg symmetries, and we shall next consider eigenvectors in another system 
of coordinates by following a more classical analytical method. 



A. Solutions through Weyl-Heisenberg symmetries 

The algebraic structure is easily displayed if we adopt the method of separation of Carte- 
sian variables. It leads to the following form of the Hamiltonian: 

* = <£ + ?"^ 2) + ( £ + \™* Y2) + T L " s H " + T L " (2 ' 3) 



where uj c = eB/mc is the cyclotron frequency, u = + AuJq, and L z = XP y — YP X . 
Here, we clearly see the splitting of the Hamiltonian into two independent harmonic 
oscillator Hamiltonians plus the angular momentum operator. Instead of directly using 
the oscillator annihilation operators: 



1 (X il „ v 1 Y il 



we work with two new ones, which are linear superposition of a x and a y 



a>d= —^(ax-ia y ), a g = ~^(a x + ia y ) , (2.5) 



where Iq — \l2%/muj. Note that a d and a g are bosonic operators: [ad, a\\ = I = [a s ,aj], 
and one has the useful identities: 



X = —{a d + a d + a g + at), Y = —{-a d + a\ + a 

L " 2l 



9 a g)i 



p x = -a\ + a g - a\), P y = -^(a d + a\-a g - a\). (2.6) 
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The operators Ti.o and L z then can be simply expressed in terms of the number operators 
N d = a\a d and N g = a g a g as: 

H = ^(iV d + iV, + l), L z = h(N d -N g ), (2.7) 

and so 

H = ^(N d + N g + 1) + ^(N d - N g ) = h(uj + N d + uj^N 9 + |), (2.8) 

where u>± — (u ± co c )/2. Eigenvalues are trivially found from the expression: 

H I n d ,n g ) = E ndng \ n d ,n g ), (2.9) 

with E ndng = h(u + n d + u_n g + uj/2) and where n d and n g are non-negative integers. The 
corresponding eigenvectors are tensor products of single Fock oscillator states: 

I n d ,n g ) =| n d )®\n g ) = -—L^aD^a^ | 0,0). (2.10) 

yjn d \n g \ 

B. su(2) and su(l, 1) symmetries 

Two distinct dynamical symmetries exist on the level of quadratic observables. 
1. su(2) symmetry 

The first one is of the type su(2) and is put into evidence by introducing the oper- 
ators 



S + = a\a g) S- = a\a* S z = -^-JL = g. (2.11) 



N d — N„ L, 

a g a d , b z = 
The commutation relations read as: 

[S + ,S-]=2S Z , [S Z ,S±] = ±S±, (2.12) 
and the invariant Casimir operator is given by 

C = \(S + S„ + S_S + ) + Si = ( ^^ )( ^i^£ + x). (2 . 13 ) 

Therefore, for a fixed value A = (n d +n g )/2 of the operator (N d +N g )/2 = Ho/(fuj) — 
1/2, there exists a (2A + l)-dimensional UIR of su(2) in which the operator S z has 
its spectral values in the range —A < c = (n d — n g )/2 < A. 

2. su(l, 1) symmetry 

The second one is of the type su(l, 1): 

1 7~t 

T + = a\a\, T_ = a g a d , T Q = -(N d + N g + 1) = ^, (2.14) 

with 

[T + ,T_] = -2T , [T ,T ± ] = ±T ± . (2.15) 
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The Casimir operator then reads as: 

V = \{T + T_+T_T + )-TZ = = (2.16) 

When n<i > n g , for a fixed value 77 = (rid — n g + l)/2 > 1/2 of the operator (N^ — 
N g + l)/2, there exists a UIR of su(l, 1) in the discrete series, in which the operator 
T has its spectral values in the infinite range 77, 77 + 1, 77 + 2, • • •. Alternatively, 
when rid < for a fixed value 1? = {—rid + n g + l)/2 > 1/2 of the operator 
(— N d + N g + l)/2, there also exists a UIR of su(l, 1) in which the spectral value of 
the operator T runs in the infinite range 1?, 1? + 1, 1? + 2, • • •. 

C. Solutions through analytical derivation 

The analytical solutions are obtained by separation of polar coordinates in the stationary 
Schrodinger equation 

W*(r, 9) = [~^(dl + -d r + - 2 d 2 e ) - i^ c d e + ^Vr 2 ]tt(r, 9) = E9(r, 9). (2.17) 
2m r r z 2 8 

We consider ty(r,9) as an eigenfunction diagonal with respect to the conserved angular 
momentum and put ty(r,0) = R(r)e tad . The function R(r) is determined in terms of 
Laguerre polynomials. Explicitly we have 

*(r,6>) = * B|a (r,6>) = (-1)" x ^L^* exp [-^] (f )H LS«»[£] e ** (2.18) 

V7T/0 V (n + |a|J! 2/5 / <o 

where n — 0, 1, 2, • • • is the principal quantum number and a = 0, ±1, ±2, • • • the angular 
moment quantum number. The eigenenergies are obtained to be: 

E na = hcu(n+ ] -^—) + -^a, (2.19) 
Therefore, n and a are related to rid and n g by the following relations: 

rid = n -\ — (M + cx), and n g = n H — (\a\ — a). 
We shall denote the eigenfunction indifferently by ^ n ,a(r, 9) — (r, 9 | n, a) = (r, | n d , n g ). 

D. Filling the shells with fermions 

Each of the above two symmetries affords a way of ordering the pairs (rid,n g ). However, 
they do not provide any hint for ordering the energy eigenvalues E nd1lg with the exception 
of two limiting cases: weak field one and strong field one. 

I. Weak field case 

Suppose the average number of electrons (N e ) obeys u c (N e ) « cu . In terms of the 
su(2) symmetry, the energy eigenvalues read as: 

E„ d „ g = hu\ + hw c q + -Hcu, (2.20) 
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and can be approximated by 



E nd n g « hcu {2\ + 1) = E\. (2.21) 

Therefore, in the weak field limit, the su(2) symmetry becomes a true symmetry 
of the Hamiltonian, which explains the degeneracy of order 2 A + 1 for the level E\. 
Note that there are (A + 1)(2A + 1) (spinless) electrons which fill the shells up to 
A . 

2. Strong field case 

In the limit of strong magnetic field oj c ^> u , we have E ndTlg « Tiu) c {n d + |). There- 
fore, for a given value of n d , we have an infinite degeneracy labelled by n g or by 
a = n d — n g < n d . The quantum number n d corresponds to the Landau level index 
(as well as n for negative a). One can reinterpret it in terms of su(l, 1) symmetry 
by noting that, for a given value of a < 0, the energy eigenstates are ladder states 
for the discrete series representation labelled by $ = —a/2 + 1/2. 

3. Generic intermediate case 

In the uncommensurate intermediate case, which means oj + /oj_ ^ Q and no approx- 
imation is relevant, we are faced to the problem of ordering the relatively dense (but 
not uniformly discrete) set of eigenenergies: 

S =^L-^L = ^± nd + n (2.22) 

In the commensurate case, u + /u- = p/q e Q, degeneracy is possible: 

E nd n g = E n i n i iff - = j. (2.23) 

d 9 q n d -n' d 



III Standard Coherent States 
A. Definitions and Properties 

The fact that the eigenvectors issued from the algebraic method are just tensor products 
of Fock harmonic oscillator eigenstates allows one to construct the corresponding coherent 
states in a standard way: 

I z ^d z n g 

| z d , z g ) =\z d )®\z g ) = exp [--(|^| 2 + \z g \ 2 )] —7=\—]= I n d' n ff) 

1 n d ,n g V n d- Jn g \ 

= exp[-^(\z d \ 2 + \z g \ 2 )]e z ^ a l | 0,0). (3.1) 

These (normalized) coherent states obey some of the usual CS properties 0: 
Pi (Eigenvector property) 

| Zdi %g) Z d | Z d , Z g ) , d g | Z d , Z g ) Z g | Z d , Z g ) . 
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P2 (Minimal uncertainty relations) 

AXAP X = ^, AY AP y = ^, 

where AX = [(z d ,z g \ X 2 \ z d ,z g ) - ((z d ,z g \ X \ z d , z g )) 2 } 1/2 , etc. 
P 3 (Temporal stability) 

e -mt/n | ZdiZg } = e -i^t/2 | e -i(w+u c )t/2 z ^ e -i(u-u c )t/2 y 

P4 (Action identity) 

H(z d ,z g ) = (z d ,z g I H I z d , z g ) = h [u + 1 z d 1 2 + cj_ I z g \ 2 + ^]. (3.2) 

The function Ti,(z d , z g ) has been called lower (resp. contravariant) symbol of the 
operator Ti by Lieb || (resp. by Berezin Q). It will plays an important role in the 
present context. 

P5 (Overlapping) 

{z' d ,z' g I z d ,z g ) = e «>(^+*^) e -(l*<i-*il 3 +l*<r-« / 9 l 2 )/2 i 
P6 (Resolution of the unity) 

1 = "9 L2 I Z d ,Zg)(Z d ,Zg I d 2 Z d d 2 Zg. 

The last property is also crucial in our context. For any observable A with suitable 
operator properties (traceclass, ...) there exists a unique upper (or covariant) symbol 
A(z d ,z g ) defined by 

A = —9 L* M z d,z g ) I z d ,z g )(z d ,z g | d 2 z d d 2 z g . (3.3) 
For instance, upper symbols for number operators are given by 

N d {z d , z g ) = \z d \ 2 - 1, N g {z d , z g ) = \z g \ 2 - 1, (3.4) 
and hence the upper symbol for our Hamiltonian ( |2.8|) 



K(zd, z g ) = Ti[u + \z d \ 2 + uj^\zg\ 2 - -]. (3.5) 
Finally, we should mention the useful trace identity for a traceclass observable A: 

TrA = — J^ 2 A(z d , z g ) d 2 z d d 2 z g = J^ 2 A(z d , z g ) d 2 z d d 2 z g , (3.6) 
where A(z d , z g ) = (z d , z g \A\ Z d i Z g ) 



B. Fock-Bargman space 

The coherent state function (r, 9 \ z d , z g ) is, up to an exponential factor, the integral 
kernel for the isometry mapping the Hilbertian span L 2 (M 2 ) of the set of eigenfunctions 
( |2 . 1 8|) onto the so-called Fock-Bargman space, i. e. the Hilbert space 



TB ^ L c 2 ntirc (C 2 , le-(l-' 2 +'^ 2 ) d 2 z d d\) (3.7) 

of entire two complex variable functions f(z d , z g ) that are square integrable with respect 
to the measure 7r~ 2 e _ (' Zd ' + ' Zfl ^ d 2 z d d 2 z g . This unitary mapping L 2 (IR 2 ) 9 ^(r, 6) — > 
Q(z d , z g ) G TB and its reciprocal are explicitely given by 

*(r,9) = \ L^r,6,z d rz 9 Mz d ,z g )e^ 2+ ^ 2 U 2 z d d 2 z g (3.8) 
$(z d ,z g ) = J R2 !C(r,-6,z d ,z g )y(r,6)rdrd6. (3.9) 
The kernel is given by the following generating function: 

rid z n g 

K{r,e,z d ,z g ) = e^ 2+ ^ 2 (r,e\z dl z g )= ^ At M ^M) 

n d ,n g Vn d l Jn g \ 



o^/Wl) e -(^-e~ W r/lo)(z g -e^r/l )_ (3.10) 



Note that normalized eigenkets ( |2.18| ) are mapped to the corresponding normalized eigen- 

states 

rid 7 ri g 

4=4^ = e^l 2+ '^/ 2 <^ g | n d ,n 9 ) (3.11) 
Vn d l Jn„\ 



in the TB representation. Also note that in this representation the operators a and a) 
have the simple form: 

d d 
a d = — , a) g = (multiplication by) z d , a g = — , a\ = z g . (3.12) 

IV Thermodynamical potential 

Let us now enter the core of the physical question we addressed in the introduction. We 
assume that the total number (N e ) of electrons is large enough for making no appreciable 
difference between a grand canonical ensemble and a canonical one. Then the magnetic 
moment M is given by 

— (I); 

where Q is the thermodynamical potential, 

fi = -^Trlog(l + e- /3(w ^ ) ), (4.2) 
P 
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with (3 = l/{ksT). The average number of electrons is given by: 



oo oo 



W = E E f( E n d n g ) = TrW) = -d,Q, (4.3) 

n,;=0 n g =0 

where f(E) = 1/(1 + e^^ E ~^') is the Fermi distribution function. The magnetic moment 
can be yielded from eqs.(fOJ) and (|4.2|) and reads as: 

M_ _2 (iV d +l/2V + -(iV g + l/2V„ 

fi B ' u 1 (i + e m-»)) 

_2 ^ ^ (n d + 1/% - (wg + 1/2 V- 
^on^o (l + /cV^ + n d +-n fl)) ' K ' } 

where k± = exp (/3(/i ± tku/2)) = k±(H,T), and /^b = he/(2mc) is the Bohr magneton. 
Despite its concise appearance, the computation of the double series (|Q ) is not eas- 
ily tractable on a numerical level. Hence we will give the preliminary estimates before 
presenting much more exploitable exact expressions. 



A. Berezin-Lieb inequalities for the thermodynamical potential 

First we note that log (1 + e - ^ 71 - - ^) is a convex function of the positive Hamiltonian TC, 
and so we can apply the Berezin-Lieb inequalities to examine the quasi-classical behaviour 
of the thermodynamical potential. These inequalities say that, for any convex function 
g(A) of the observable A, we have 

J c2 g(A)d 2 z d d 2 z g < Tig(A) < ± g(A) d 2 z d d 2 z g . (4.5) 

Applying them to the (concave) thermodynamical potential leads to the inequalities: 

" b le log (1 + e ' m '^ d " Zd d2 ^ - n -~]bIc logil + e ~ m ~^ dhd d2z *- 

(4.6) 

Let us insert eqs.( ^.2[ ) and Q3.5| ) and perform the angular integrations. It leads to the 
inequalities: 

1 /"OO roo 

- - / du d / du q log (i + e -/^(™+"-%-f )-*«)) < n 

p Jo Jo 

1 roo roo 

n<-o/ i g( 1 + e -/3w- + « d +--« 9 +f)-^)) ; (4.7) 

p Jo Jo 

where u d = \z d \ 2 and u g = \z g \ 2 . Changing the integration variables, say u = j3h(u + u d + 
uj_u g ), v = /3hu + u d , performing an integration by part, and introducing the control pa- 
rameters k± defined in (|4.4p , we easily reduce the above inequalities to the following 
ones: 

<P(k+) <tt< (4.8) 
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where the function <j> is given by: 

K 



oo u 2 e ~ u 



2(3(phu ) 2 Jo 1 + ne~ u 
I 

(log «) 3 



du 



/3(/3tLo) 2 F ^ ^ 
1 

/3(/3^o) 2 



7T log K 



for k < 1, 
for ft > 1. 



(4.9) 



We have introduced here the function F s , of the Riemann-Fermi-Dirac type, defined as: 



Fs(z) 



oo 



E 

m=l 



(4.10) 



In the high temperature region |// ± Ku/2\ <C /cgT we have k± « 1. From (|4.8|) we see 
that the thermodynamical potential is approximately equal to: 



F s (-1) « -0.901543 A; B T 



We now consider the more realistic case: /i ^> /k<j/2 and [i ^> k B T. Let us split the 
function <p into three parts: 

(f)(K ± ) = At^- + s±, 



(4.11) 



with 



A = 


n 

2 


A 








s± = 


k B T 



1 / // 



3 V^cuo 
1 / /i ^ 2 



2 1/cj\ 2 7l 2 k B T" 
' ' ' " ' IT '0 



4 VuV 



2 V/icuc 
/ k B T\ 



+ 



u; Y 



7T 2 /fc B T 



24 Vcuq/ 6 \ ftcuo 



F 3 (-exp [-/%±^/2)]). 



(4.12) 



Then we see that Q lies in the interval [A + S + — A/2, A + S_ + A/2]. Replacing S± by 
the approximate expression 



5*0 = k B T 



k B r 



(4.13) 



and observing that the ratio 



A 



Too 



\A + Sb| /i + 7r 2 fkjjzy + | _ (hMLY F ^_ e -p»\ 
tends to zero, we see that the thermodynamical potential can be estimated as: 

n « a + s 



(4.14) 



2 



3 \TlujqJ 4 Vcuq/ 3 \ nu Q 



+ k B T 



k B r 



F 3 (-e-^)-(4.15) 
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A similar asymptotic behaviour holds for the thermodynamical limit (N e ) — ► oo. Indeed, 
in this quasiclassical regime, the average number of electrons is calculated to be: 



(N e ) « -dJA + S t 



o 



1 ( I 1 



> 

for /i > k B T and /i > hw/2. (4.16) 



2 V^cuo 

The magnetic moment hence reads as: 

M = Xp H, with Xp = V ■ ( 4 - 17 ) 

Hence, we can assert that the system shows an orbital paramagnetism in this limit. At 
this point, we refer to a recent work by Combescure and Robert || in which precise infor- 
mations are given for the magnetisation of electron gas constrained by general confinement 
potentials. 

B. Exact expressions for the thermodynamical potential 

We now determine the thermodynamical potential in a precise way by applying the for- 
mulas ( |7.5| ) and ( [7.7| ) in the appendix. The function Q(k) defined by ( |7.6j ) takes the 
following closed form: 

e(Jfe) = Tr(e-( lfc+1 ^ w ) = e ~ (ik+1) ^ 1 — B 1 — a . (4.18) 

The Fourier integral representation for the thermodynamical potential hence reads as: 

pJ-oc 2coshffc \ik + l) Vi _ e -(ifc+i)f y \i _ e -(tfc+i)f»«-J 

(4.19) 

As indicated in the formula ( |7.6|) , this Fourier integral is given as a series by using the 
residue theorem. One can easily see that the numbers (2m + m G Z are simple poles 
of sech^fc, i is a double pole of 0(fc), and i + 47rm/ (j3hu + ), % + 4irm/ (/3hu)-), m e Z* are 
simple or double poles of 0(/c) according to whether u + and u;_ are uncommensurable 
or not (see Fig.l). In order to fulfill the requirements of the Jordan Lemma, one has to 
consider the following two cases: fi < huj/2 and \x > Tiuj/2. In the first case we take an 
integration path lying in the lower half-plane and involving only the simple poles (2m+l)z, 
m < 0. It leads to the result: 



tl = ±ST ^ 3 , , (4.20) 

4p m=1 m sinh (|^,o; + m) sinh (^huj-m) 



which corresponds to the case k < 1 in Eq.( [4.9|) . In the second case, an integration path 
in the upper half-plane is chosen. It encircles all the other poles: (2m + m > 0, 
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i + 47Tm/(l3hu) + ), i + 4irm/(/3hu)-), m G Z*, as shown in Fig.l. We present the result in 
a manner which will render apparent the various regimes: 



n 




+ 



02 



+ 



^osc 



' 1 4tt 

a-i((2m + 1)2) + £ (a_!(z + — — 

m>l m±^0 " ^ 



The first term is at the origin of the Landau diamagnetism: 

2 1 



UJ, 



JL 

24 \u Q 



m±)). 

(4.21) 

(4.22) 



where xl 



|/i (^^J = — \Dqii 2 b is the Landau diamagnetic susceptibility. The coeffi- 
cient D = fi/(hu ) 2 can be interpreted as the density of states at Fermi energy. Notice 
that the value of xl is equal to one third of the "quasiclassical" one Xp obtained in ( |4.17| ), 
a feature which is reminiscent of which we encounter in the studies of free 3D electron 
gas. The second term, which gives no contribution to the magnetization, is written as: 

2 



01 



6 



( » 



\Tiojc 



k B T 



The third term is given as: 

^02 = 



J_ ~ (~1) ? 



exp 



k B T 



m) 



sinh 



m) sinh 



-m) 



(4.23) 



(4.24) 



2k B T " v 2fcflT 

It becomes negligible at low temperature regime k B T <C \i. The sum of Ql and Qqi 
is analogue to the term A in Eq.( |4.11| ) and Qq2 corresponds to S±. The last term is 
responsible for the oscillatory behaviour. We need to distinguish between irrational values 
of u) + /lo_ and rational ones: 



case oj + /oj- £ 



a 



2/5 tx 



in 



sin (^-7rm) 



sin (— Tim 



sinh (^ E ^7r 2 m) 



+ 



sin (j|jf-7rJ7i) 



-7rm) sinh ( 2 ^ BT TT 2 m) 



sin 



case lu + /l}- = p/q G 



1 

2^ 



E 

m=l, m^O mod q 



gcd (p, g) 
f-l) m 



(4.25) 



1, u + /p = uu_/q = 21 /{hP) G 



sm^vrmj 



sin(— 7rm) sinh ( 2 ^ T Tc 2 m) 



m 



+ 



+ 



E • 

m=l, m^O modp 



sin i 



Tim) 



m sin (^7rm) sinh (^^7r 2 m) 



(4.26) 



1 



E 



(_l)(p+g)fc 



ipg t^x /csinh(^-fc) 



/i /i7rA; 
cos I 



k B Tl 



/ , z 71 " 2 , x I \ ■ , ^k , 
( 7 rcoth(-A;) + -)sin(— 
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C. The thermodynamical potential via symbol calculus 

Let us define the thermodynamical potential operator as O = — -s log (1 + exp [— (3(Jht — /i)])- 
We have: 



° = ~2 L 0( z d,Zg) I z d ,z g )(z d ,z g I d 2 z d d 2 z. 



0{Zd,Zg) — (Z d ,Zg l O l Z d ,Zg), 



(4.27) 



We then perform the angular integrations and take its trace. We get an integral repre- 
sentation of the thermodynamical potential Q = TrO: 



oo roo 



| roo roo A j 

Q = -- / / du d du g O(u dj Ug) = -— . , 
p Jo Jo p Jo Jo 



du d du g 0(u d , u g ). 



(4.28) 



The problem turns out to evaluate the upper and lower symbols. The lower one is shown 
here with the integral representation: 



6{u d ,u 



(u d +u g ) r+OQ _-(iAH-l)£(^-/i) 



2(3 



(cosh %k){ik + 1 



exp 



(u d e~ {lk+l) 2 hw+ ) exp ( M9 e- (tfc+1) i^-) rfA; 



(4.29) 



V Average number of electrons and magnetic moment 



In this section, we will exploit the formulas ( |4.20| )-( |4.26| ) to obtain the exact expressions 
of the average number of electrons and the magnetization. We will restrict ourselves to 
the more realistic case: // < Too/ 2. 

The average number of electrons is easily derived by taking the derivative of —Q with 
respect to \i. It is found to be: 



(Ne) 



24 W / 

1 oo 



+ 



+ 



7T 



m=l 

oo 



sinh (^tuu + m) sinh (| ftiO- 



m) 



m=l 



cos (-j^-irm 



hu_ sin (^vrm) sinh (f^vr 2 



+ 



cos (-^jj-nm) 



m) 



%uj + sin (^vrm) sinh (^^ir 2 m) 



= (N e ) L + (N e ) 01 + (N e ) 02 + (N e )- c + (N e ) 



+ 

osc" 



(5.1) 



The magnetic moment is decomposed into four parts and is expressed in Bohr mag- 
neton units: 



M 



02 



dhuj r 



2/i 



B 



dhidr 



= 2u jB {M l + M + M osc + M+ c ), 



(5.2) 
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where 

L Y2%uq \ojq) 2fi B ^ L ' ^ ^ 

1 ~ m [g; + coth (phuj + m/ 2) - u„ coth (phuj_m/ 2)] 

~ J sinh(/3^ + m/2)sinh(^_m/2) ' 1 ' 

and, for the irrational case uj + /oj_ ^ Q, 

M - = ksT f (-irsin(27rm/i/(^_)) y 

Tiuj sin (irmLJ + /u-) sinh {2-K 2 mk B T / {TiujJ)) 

7TU / |i \ 7TCo> + / CJ + \ 7T 2 k B T ( 2 /c^T' 

cot 2i\m- cot 7rm — coth 2ix m- 



M+ = k B T « (-l) m sin(27rm/i/(^ + )) ^ 

Ttu; , ( ~ L sin (7rmct;_/c(j + ) sinh (27r 2 m/cB^/(^+)) 

7TU / u \ 7ru;_ / c<j_\ n 2 k B T ( 2 fc#T' 
cot znm- cot nm — coth 2n m- 



We will not give the expressions of J^A^sc * n the rational case because the magnetization is 
a continuous function of uo c and its behavior can be fully understood from the irrational 
one. 

Discussion 

The temperature scale is compared to the two natural modes u± of the system and draws 
three possible intrinsic regimes: high temperature regime k B T > Huj + , low temperature 
regime k B T < huj_, and intermediate temperature regime huj_ < k B T < hu> + . Remember 
that we work in the large electron number region: /i > Tiuj/2. 

1. High temperature regime: k B T > Huj + > Tiuj-. 

This inequality implies the following constraint on the field: 

Wc < k B T _ hu ^ k B T ^ ^ 

Uo hujQ k B T Tiujq 

We can see that M.q is the dominant term for the magnetic moment in regard to 
M. osc because of the smallness of arguments of the sinh (in the denominator) and 
coth (in the numerator) functions. Hence, M pa 2/j b (M.l + -Mo), which shows 
mainly Landau diamagnetism. Similarly, we infer from (|5.1|) that (N e ) ~ (-We)z, + 
(N e ) 01 + (N e ) 02 . 

2. Low temperature regime: k B T < hu^. 

The magnetic field is restricted by the inequality: 

< tiuo_ _ ksT _ hwo_ ^ ^ 

ujq k B T huj k B T 



14 



Now the M.q term becomes excessively small due to the rapidly decreasing expo- 
nential factor and the large arguments of the hyperbolic functions present in the 
expression. The magnetization is hence approximately determined by the three 
terms: Ml, Mtsci an d M~ sc and exhibits oscillating behavior. 

(a) Strong fields uo c 3> uoq. 

We make the following approximations: 

w+ »,Ji + (i*YY w _»5£, ^(^) 2 . (5.9) 



^«i±^+lf^) . (5.10) 



Assume also that uo c < 2u y (fi/h) 2 — 1 to insure the validity of /i > huj/2. One 
can see that the denominator of -M^ sc contains the product of sin and two sinh's 
with small argument and, hence, M ~ 2jj, B {M. L + M-tsc)- I* 1 particular, after 
replacing u + by uj c in the sinus argument of the numerator of -Mq SC , we find 
that the magnetization is periodic with respect to the inverse of the magnetic 
field, a characteristic fact of the "de Haas-van Alphen" regime. The similar 
behavior holds for the electron number: (N e ) ps (N e ) l + (fJ, / huo) 2 / 2 + (N e )+ c . 

(b) Weak fields uj c <C cu . 

In this case the two characteristic frequencies and their ratio are approximated 
to 

, , 1 UJ C 1 / co c \ 2 

1 ± — - + - — 

2uq 8 \ujqJ uj t ujq 2 \ujo< 

Now M~ sc and M.„ sc have the same order of contribution due to the presence 
in their denominators of 

(sinh (2tx 2 mk B T /huJ T )f « (sinh (27r 2 mk B T/hu )) 2 (5.11) 

and of the strongly oscillating functions 

(sin {jmuj±/uJ T )) 2 ~ (sin (7imu c /iJ )) 2 . (5-12) 

Consequently, the system is considered as lying in the mesoscopic phase. Simi- 
lar conclusions can be reached for the behaviour of the average electron number. 

3. Intermediate temperatures: hu^ < k B T < huj + . 
These inequalities imply the constraint 



UJ C 

— > 

u 



TlUJa knT 



k B T huo 



(5.13) 



The weak field case occurs only when k B T approaches hu . In this case, one may 
think that the oscillatory terms JA~ SC and /A„ sc give their contribution to M as 
what we have seen in the previous subsection. But, in fact, the hyperbolic sinus 
functions of the denominators have large arguments: 

sinh (2n 2 mk B T /hcu T ) ^ sinh(19.74m) ^ 0.5 exp(19.74m), 



15 



and so overcomes the algebraic contribution of the sinus functions. The system 
hence goes to the Landau diamagnetic regime. On the other hand, for a strong 
field, we return to the approximation M « 2[i b (Ml + A4 + J and the system shows 
the de Haas- van Alphen effect. 

We have repeated in Fig. 2 the phase diagram of magnetization proposed by Ishikawa and 
Fukuyama. Let us however mention that the two intrinsic frequencies u± do not represent 
the exact borders between the different magnetic phases. It should be taken in a qualita- 
tive sense only. In order to justify this, we choose the chemical potential equal to 100.0h.ujQ 
and make temperature vary through the different magnetic field regimes. Fig. 3 (a) illus- 
trates the weak field regime. One can see that at low temperature ksT = 0.00lhui the 
magnetization experiences strong fluctuations, called mesoscopic fluctuations. As T in- 
creases (ksT = O.lhuJo), the strength of these fluctuations decreases. As the temperature 
is getting higher (for example for ksT = 0.5hujo), the fluctuations disappear and we reach 
the Landau diamagnetism. In Fig.3(b) the magnetic field lies between l.^uj^mc/e and 
3.1uj mc/e. When ksT = 0.01Tiuj , the system shows large fluctuations. One can see that, 
if the magnetic filed increases, the fluctuations diminish and the cycloid-like curve appears 
(de Haas-van Alphen oscillations). At higher temperature, e.g. ksT = O.lhuo, one can 
see more clearly the phase change from the mesoscopic fluctuations to the de Haas- van 
Alphen oscillations. If the temperature continues increasing (for example, k B T = 0.5hu , 
k B T = l.Otuvo), the system will go from the Landau diamagnetism to the de Haas-van 
Alphen oscillations. Fig.3(c) shows the de Haas-van Alphen phases. The position of the 
peak can be predicted from the simple formula: 

^ c 2u nhuo . . . 

— = — , where n is some odd positive integer. 

ojq nnojQ 2 ii 

In Fig. 3(d) we see that, at the extreme low temperature {ksT = 0.001hu ), the curve 
shows some "width", which comes from the limit of the picture resolution and which 
corresponds, in fact, to small fluctuations. This gives us an example of de Haas- van 
Alphen oscillations mixed with mesoscopic fluctuations. 

VI Current Distribution 

We now turn our attention to the spatial density of current. It can be obtained by the 
following formula: 

J(r) = ft ^t( R ) tfl (p + ^A(R)) ^(R)^ , (6.1) 

where (•) = Tr(/(7Y)(-)) is the thermal average, and ip(R) the field operator: 

^(R) | n,a) = * n ,a(r). 
Due to the symmetry of the system, the current distribution is purely orthoradial, 

J(r) = J e (r)e e , (6.2) 
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and its component is given by the series: 

£ w e r 



J e {r) = -ev £ (c4 + f ^(r, 9)f{E na ), (6.3) 



where £ = yh/(mu ) is the characteristic length and = uj £ the characteristic speed in 
the harmonic potential. The current distribution induces a local magnetic moment, 



M(r) = — r x J = j- r j e (r)e z = M 2 (r)e 2 



which in turn gives rise to the total magnetic moment, 



M 



M z {r) dS=- 

c Jo 



J@(r)r 2 dr. 



By applying Eq. (|3.8| ) to the expression (|6.3j ) of Je(r) we get 



(6.4) 



(6.5) 



Mr) 



- e ) v o f(Hu + n d + u„n g + — ) I (n d 



n. 



n d ,n g 



uj c r 



n' 



to)"- ( 



d 2 z 1 ^ d^Zg d 2 z d d 2 Zg. 



(6.6) 



The above equation is rotational invariant and one can hence put 6 = 0. Let us now 
represent the Fermi-Dirac function by the Fourier integral ( |7.2| ) and perform the discrete 
summation. The four integrals on the complex plane can be easily performed by using 
well-known gaussian integrals of the type 



°° e - t2+tz dt = ^e^(-Z 2 ). 

-oo 4 



We finally obtain the following expression: 



Mr) 



-evo r 



00 dk exp [(ik + 1) § //] 

4 > 



16ttZ§ £ 7-oo cosh (f fc) sinh 2 ((iA; + l)/3 

^ sinh ((iJfe + 1) A - — sinh ((ik + l)/3^) 
cjq 4 cuq 4 ' 

-2r 2 sinh((i£; + l)/3^) sinh((ik + l)/3^) 



x 



exp 



11 sinh ((ik + l)/3 



(6.7) 



Three types of singular points appear in the integrand: 

1. Pole k = i. 

This pole is triple since it appears once in cosh and twice in (sinh (ik + l)/3 

2. Poles k = (2m + m e Z*. 

They are the other (simple) poles of cosh |/c. 
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3. Singular points k = i + 47rm/ (ftfuv), m G Z*. 

These points are essential singularities since they appear in the exponential. 

The current distribution is hence separated into three components responsible for the 
Landau diamagnetism, weakly diamagnetism, and oscillation behavior: 



In a forthcoming paper, we shall go deeper into the analysis of such an expression for the 
current distribution. 

VII Conclusion 

In this paper, we have established exact formulas for the thermodynamical potential of a 
two-dimensional gas of spinless electrons, confined in an isotropic harmonic potential, and 
submitted to a constant perpendicular magnetic field. From this it has become possible 
to study the magnetic moment and other thermodynamical quantities at different regimes 
of temperature and field. This exhaustive study was made possible thanks to the specific 
simplicity of the isotropic harmonic potential. Of course, there exist other situations in 
which we could get similar exact expressions: anisotropic harmonic potential, harmonic 
groove, as indicated by 0. In a next future, we shall deal with less tractable but 
still integrable models (see for instance 0, like infinite cylinder potential or quantum 
rings, cases in which the expressing of the trace (7.6) in a closed form is quitely not 
expected. We shall also intend to make use of other families of coherent states, possibly 
more adapted to these new situations, in order to use an efficient symbol calculus. For 
instance, for the infinite rectangular well potential, we are thinking about using tensor 
products of infinite well coherent states introduced and analysed in a recent paper 
Finally, we shall also explore the behaviour, at different temperature and field regimes, of 
other thermodynamical quantities of current experimental interest, like the heat capacity 



Appendix: Fermi-Dirac trace formulas 

It is well known that, like the Gaussian function, the function sechx = 1/ coshx is a fixed 
point for the Fourier transform in the Schwartz space: 



Hence, given an Hamiltonian H,, we can write for the corresponding Fermi operator: 




PH. 




(7.1) 



1 



(ifc+l)f(H- M ) 



W) 



1 + e/S(«-M) 



4 cosh | k 



dk. 



(7.2) 



— oo 
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Similarly, we can write for the thermodynamical potential operator: 

1 a ,~. x 1 r+°° e -(^+i)f (w-m) 

" p l0E I 1 + e ~ m ^ = -Jj L (2co 3h p)( rt+ l) *" < 7 - 3 ' 

Therefore, the average number of fermions and the thermodynamical potential can be 
written (at least formally) as follows: 

/+oo 

1 1 r+oo P (ik+1)@£ 

where designates the function 

e(Jfe) = Tr(e" (ifc+1) ^). (7.6) 

Observe that (2m + m e Z are (simple) poles for the function 1/ cosh |& and z is 
a pole for the functions Q(k) and 1 / (zA; + 1). These Fourier integrals can be evaluated 
by using residue theorems if the integrand functions $i(&) = Q{k)/ cosh|/c and $2(&) — 
e(k)/((ik + 1) cosh | A;) satisfy the Jordan Lemma, that is, § x {Re i6 ) < g(R), § 2 {Re i6 ) < 
h(R), for all 9 G [0, 7r], and g (-R) and vanish as -R — > oo. The quantities (N) and Q 
are then formally given by 



2ni 



(7.7) 



+ E i-i((2m + l)») + E 11 -iW 

m=l v 

where a-i(-) denotes the residue of the involved integrand at pole (■), and the fc„'s are 
the poles (with the exclusion of the pole i) of Q(k) in the complex fc-plane. 

We now introduce the spectral resolution of the (bounded below) self-adjoint operator 

H: 

r+oo 

<p(H)= <p(X)E(dX), (7.8) 



where p is a complex- valued function and E\ = JL^ E(dX) is the resolution of the identity 
for the Hamiltonian Ji. Define the density of states z/(A) as TrE(dX)/d\. The trace 
formula ensues: 

r+oo 

TMH) = / <p(X) u(X) dX. (7.9) 



Let us now introduce the weighted density of states w(X) = e 2 X u(X) and its Fourier 
transform 

j r+oo 

w{k)^-f= e~ lkX w(X)dX. (7.10) 



/2vr 

Then, from (|7.4j ), ( [7.5|) and (|7.9|) , we can represent (N) and fl as follows: 



, r+oo e (ik+l)^ a ^ 



/3 



2tt e( ife+1 ^ A/ /3 n J( 2tt 



0u 



a = -irL w^m^f ^ k)dk= -T> e (7 ' 12) 
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where we have introduced the weighted functions Z\ = seeh(^k)w(k) and Z 2 = sech(7r/c/ (3){i2k/ (3+ 

iy l w(k). 
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Figure caption 

Fig.l Poles of the Fourier representation of the thermodynamical potential Q. The poles 
lying on the imaginary axis are simple except for the point % which is of order 4. 
The poles lying on the line k = i in the complex A;-plane (i is excluded) may be 
simple or double depending on whether u + and are uncommensurable or not. 

Fig. 2 Phase diagram of the magnetization. In high temperature and low magnetic field 
region, the system shows Landau diamagnetism; in low temperature and low mag- 
netic field region, mesoscopic fluctuations appear; and in strong magnetic field re- 
gion, the system experiences the de Haas-van Alphen oscillation phase. The two 
curves, k B T = hu + and k B T = hu-, give a qualitative indication about the phase 
borders. 

Fig. 3 Magnetization curves versus the magnetic field at different temperatures. We rep- 
resent the cyclotron frequency in uj q units and the magnetization in 2//^ units. The 
chemical potential \i is set up to 100.0fru;o. 

(a) oj c is less than 0.8u Q . Temperatures are chosen to be 0.001hu , 0.1hu , and 
0.5hu , respectively. 

(b) oj c is between 1.9uj q and 3.1c<j . Temperatures are chosen to be O.Ol^o, 0.1^o;o, 
0.5hujQ, and 1.0huj , respectively. 

(c) oj c is between 4.0cc;o and 15.0u;o- Temperatures are chosen to be O.l/iwo, 0.5hu , 
I.OHcuq, and B.OHcuo, respectively. 

(d) uj c is greater than 15.0o;o- Temperatures are chosen to be O.OOl/k^o, 0.5/iu;o, 
l.Ohcuo, and 5.0hu , respectively. 
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